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We present an observation of the decay B° — > -k°-k° based on a sample of 124 million BB pairs 
recorded by the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We observe 
46 ± 13 ± 3 events, where the first error is statistical and the second is systematic, corresponding 
to a significance of 4.2 standard deviations including systematic uncertainties. We measure the 
branching fraction B(B° -> ttV) = (2.1 ± 0.6 ± 0.3) x 10~ 6 , averaged over B° and B° decays. 
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PACS numbers: 13.25.Hw, 11.30.Er 12.15.Hh 

The study of B meson decays into charmless hadronic 
final states plays an important role in the understanding 
of CP violation in the B system. In the Standard Model, 
CP violation arises from a single complex phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix V lj. Measurements of the time-dependent CP- 
violating asymmetry in the B — > n + ir~ decay mode 
by the BABAR and Belle collaborations provide in- 
formation on the angle a = arg [— V td V r t ^/V ud l / u * b ] of the 
unitarity triangle. However, in contrast to the theoret- 
ically clean determination of the angle (3 in B° decays 
to charmonium plus neutral- kaon final states 0, 0] , the 
extraction of a in B° — > tt + it~ is complicated by the 
interference of amplitudes with different weak phases. 
The difference between a e g, derived from the measured 
B a — > 7r + 7r - asymmetry, and a may be evaluated using 
isospin relations between the amplitudes for the decays 
B °(B°) 7T+7T-, B°(B°) -> 7r°7r°, and B ± -> t^tt @. 

The primary contributions to the decay B° — > 7r 7r° 
are expected to come from the so-called color-suppressed 
tree and gluonic penguin amplitudes |(J. The branching 
fraction for B° — > 7r°7r° has been calculated in various 
QCD models 0. All models use as inputs the values of 
the CKM angles, typically taken from unitarity-triangle 
fits. The predictions for B(B° — > 7r°7r°) are in the range 
(0.3-1.1) x 1CT 6 . 

In this paper, we report the observation of the decay 
B° -> 7r°7r° based on (124 ± 1) x 10 6 T(45) -> BB pairs 
(on-resonance) , collected with the BABAR detector. We 
also use approximately 12 fb -1 of data recorded 40MeV 
below the BB threshold (off-resonance). 

BABAR is a solenoidal detector optimized for the 
asymmetric-energy beams at PEP-II and is described 
in detail in Ref. Charged particle (track) momenta 
are measured with a 5-layer double-sided silicon vertex 
tracker and a 40- layer drift chamber inside a 1.5-T su- 
perconducting solenoidal magnet. Neutral cluster (pho- 
ton) positions and energies are measured with an electro- 
magnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals. The photon energy resolution is ue/E = 
|2.3/ J B(GeV) 1 / 4 © 1.9} %, and the angular resolution 
from the interaction point is gq = 3.9°/y /, E( GeV). 
The photon energy scale is determined using symmetric 
7T° — * 77 decays. Charged hadrons are identified with 
a detector of internally reflected Cherenkov light and 
ionization in the tracking detectors. The instrumented 
magnetic-flux return detects neutral hadrons and iden- 
tifies muons. High efficiency for recording BB events in 
which one B decays with low multiplicity is achieved with 
a two-level trigger with complementary tracking-based 
and calorimetry-based trigger decisions. 

Candidate n° mesons are reconstructed as pairs of pho- 
tons, spatially separated in the EMC, with an invariant 



mass within 3tr of the ir° mass. The mass resolution a is 
approximately 8 MeV/c 2 for high-momentum tt° mesons. 
Photon candidates are required to be consistent with the 
expected lateral shower shape, not to be matched to a 
track, and to have a minimum energy of 30 MeV. To re- 
duce the background from false n° candidates, the angle 
9-y between the photon momentum vector in the tt° rest 
frame and the ir° momentum vector in the laboratory 
frame is required to satisfy | cos# 7 | < 0.95. 

B meson candidates are reconstructed by combin- 
ing two 7T° candidates. Two kinematic variables, used 
to isolate the B° — * 7r°7r° signal, take advantage of 
the kinematic constraints of B mesons produced at the 
T(4S). The first is the beam-energy-substituted mass 
toes = \/( s / 2 + Pi ' PB) 2 /Ef - p|, where yfs is the to- 
tal e + e~ center-of-mass (CM) energy. (Ei, Pi) is the four- 
momentum of the initial e + e~ system and ps is the B 
candidate momentum, both measured in the laboratory 
frame. The second variable is AE = Eb — y/s/2, where 
Eb is the B candidate energy in the CM frame. The AE 
resolution for signal is approximately 80 MeV. 

The primary source of background is e + e~ — > qq (q = 
u,d,s,c) events where a ir° from each quark jet randomly 
combine to mimic a B decay. The jet-like qq background 
is suppressed by requiring that the angle 9 S between the 
sphericity Q axes of the B candidate and of the remain- 
ing tracks and photons in the event, in the CM frame, 
satisfies | cos# s | < 0.7. The other source of background 
is B ± — ► p ± 7r° ( / o ± — * ^n ) decays in which the charged 
pion is emitted nearly at rest in the B rest frame so that 
the remaining two tt° mesons are kinematically consistent 
with a B° — > 7r°7r° decay. Energy resolution smearing 
causes some B ± — > p ±r K Q events to have AE above the 
kinematic limit of tub — m^. From simulation, other B 
decays contribute no more than one background event. 

The number of signal B° — > 7r°7r° candidates is deter- 
mined in an extended, unbinned maximum-likelihood fit. 
The variables used in the fit are toes, AE, and a Fisher 
discriminant F. The F discriminant is a linear combi- 
nation of three variables, optimized to separate signal 
from qq background. The first two variables are sums: 
Lq = J2iPi an d L2 — J2iPi cos2 @i where pi is the mo- 
mentum and 0i is the angle with respect to the thrust 
axis of the B candidate, both in the CM frame, for all 
tracks and neutral clusters not used to reconstruct the 
B meson. The third variable in F is the output of a 
neural network designed to separate B events from qq 
background, whose inputs are information from the re- 
maining tracks and photons in the event. The inputs 
include information about high- momentum leptons, low- 
momentum leptons, charged kaons, and slow pions (from 
D* + — > D °7r s | ow ) in the event; these are the same inputs 
used in the B-tagging algorithm of Ref. Q. All neural- 
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network training and Fisher-discriminant optimization is 
performed using simulated events. 

The data are divided into two samples: a signal sam- 
ple with candidates satisfying toes > 5.2GeV/c 2 and 
\AE\ < 0.2 GeV, and a sideband sample with candi- 
dates from on-resonance data with toes > 5.2 GeV/ c 2 
and 0.2 < \AE\ < 0.4 GeV (and well outside the triangu- 
lar region in toes & n d AE populated by B ± — > p^n 
decays) and candidates from off-resonance data with 
toes > 5.2GeV/c 2 and \AE\ < 0.4 GeV. The sideband 
sample contains only qq background candidates and is 
used in the fit to improve the statistical precision of the 
F distribution for qq events. There are 4470 events in 
the signal sample and 3253 events in the sideband sam- 
ple. The reconstruction efficiency for B° — > ir°ir° is 
(17.7±2.7)%, and for B ± -► p ± n° is (0.8±0.1)%, derived 
from simulation. The errors are due to a systematic un- 
certainty in the efficiency for high-momentum ir° mesons 
to pass the selection criteria. 

For candidates in the signal sample the probabilities 
Vi(xj]di) used in the maximum- likelihood fit are the 
product of probability density functions (PDFs) for the 
variables Xj = {toes, AE, F}, given the set of param- 
eters oii. The likelihood function is given by a prod- 
uct over all j = 1 — N candidates and a sum over the 
i = {B° -> 7r°7r°, B ± -> p ± ir°, qq] hypotheses: 



JV 



£ = exp -$^n< [] 



^riiVi (xf,cti) 



(1) 



The coefficients n; L are the numbers of B° 



7r°7T sig- 



nal, B d 



: 7T° background, and qq background events 



p ± n° events is 



in the sample. The number of B 
fixed in the fit to the expected value based on the mea- 
sured branching fraction B(B ± — > p ± -K°) = (11.0±2.7) x 
10~ 6 10]. For candidates in the sideband sample, the 
likelihood function includes only the PDF for the F vari- 
able, and only the component for qq background. A si- 
multaneous fit to both signal sample and sideband sample 
data is performed. Monte Carlo simulations are used to 
verify that the likelihood fit is unbiased. 

The PDFs are determined from data and simulation. 
The toes an d AE variables are correlated for both B° — > 
7r°7r° and S ± — > /9 ± 7r°, so a two-dimensional PDF derived 
from a smoothed, simulated distribution is used. The 
toes distribution for qq events is modeled as a thresh- 
old function 0] whose shape parameter is determined 
from data with | cos# s | > 0.9. The AE distribution for 
qq events is modeled as a quadratic polynomial with pa- 
rameters determined from data with toes < 5.26 GeV/c 2 . 

The PDF for the F variable is modeled as a paramet- 
ric step function (PSF) for B° -> 7r°7r°, B ± -> p ± ir°, 
and qq events. A PSF is a binned distribution (as in 
a histogram), whose parameters are the heights of each 
bin. Since the parent distribution of F is not known, any 
functional form (such as a multiple Gaussian) assumed 



for the PDF will suffer from a systematic uncertainty 
due to the choice of function. By binning the data, the 
PSF substantially reduces this uncertainty. The PSF 
is normalized to one, so that the number of free pa- 
rameters is the number of bins minus one. For both 
B° -> 7r°7r° and B ± -> P ± tt°, the F PSF parameters 
are taken from a sample of 3.2 x 10 4 fully reconstructed 
B° -> L>(*W(n = 1,2,3) events in data. The F PSF 
has ten bins, with bin limits chosen so that each bin con- 
tains approximately 10% of the B° — > D^'nir events. 
Simulation is used to verify that the same distribution 
can be used for both B° — > 7r°7r° and B ± — > p^n . For 
qq background, the F PSF parameters are free parame- 
ters in the fit; these parameters are determined from data 
in both the signal and sideband samples. 

All event-selection requirements, PDF parameters, and 
maximum-likelihood fit conditions were determined be- 
fore fitting the data. 

The result of the fit is ubo^^o^q = 46 ± 13 events, cor- 
responding to a branching fraction of B(B° — > 7r°7r°) = 
(2.1±0.6) x 10~ 6 . B° and B° decays are not separated, so 
the branching fraction is measured for the average of B° 
and B°. The toes, AE, and F distributions are shown 
in Fig. ^ f° r a U data used in the fit, and in Fig. [21 for 
events that pass a requirement on the signal probability 
ratio. This requirement is optimized to maximize the ra- 
tio S/ y/S + B, where S is the number of signal events 
and B is the number of background events in the plot. 
The significance of the event yield is evaluated from the 
square root of the change in —2 In £ between the nominal 
fit and a separate fit in which the signal yield is fixed to 
zero, and is found to be 4.7ct with statistical errors only. 

The number of signal events is stable when the qq Toes 
and AE PDF parameters, or n B ±^ p ±^o , are allowed to 
vary in the fit. A validation of the maximum- likelihood 
fit is made by performing a simpler event-counting analy- 
sis, based on the number of events satisfying tighter Toes, 
AE, and F requirements. The event-counting analysis 
finds 13 ± 6 events with an efficiency of 31% relative to 
the maximum-likelihood fit. This agrees well with the 
fitted result, and has a statistical significance of 2.7er. 

This result is consistent with our previous limit for 
this decay based on 88 x 10 6 BB pairs. The data de- 
scribed in Ref. ^3] were reanalyzed with improved EMC 
energy calibration and tracking alignment. More events 
are observed in this data sample after the reanalysis, con- 
sistent with the improved understanding of the detector. 

Systematic uncertainties on the event yield are eval- 
uated by varying the fixed parameters and refitting the 
data, and are summarized in Tabled The shape param- 
eter for the threshold function describing the toes distri- 
bution for qq events is varied to account for the statistical 
error from the fit to the sample with | cos# s | > 0.9 and 
the extrapolation from | cos(9 s | > 0.9 to | cos# s | < 0.7. 
The qq AE polynomial parameters are varied by their 
statistical errors. The number of B^ — > p^n back- 
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ground events is varied according to the uncertainties 
on the B^ — > branching fraction and reconstruc- 

tion efficiency. Finally, the uncertainty in the mean of 
the AE distribution for B° — > 7r°7r° is evaluated from a 
study of B ± — ► /^tt events that have a high momentum 
7T°. Extrapolating from the uncertainty in the mean AE 
for this sample, we vary the mean of AE by ±12 MeV to 
evaluate the systematic uncertainty on the signal yield. 
The effect of these uncertainties on the significance of the 
event yield is evaluated by choosing the variation that re- 
duces the signal in all four systematic effects, and then 
refitting the data. The significance of the signal yield af- 
ter accounting for systematic uncertainties is 4.2(7. The 
change in — 21n£ as a function of the signal event yield 
is shown in Fig. 0I. 

We observe 46 ± 13 ± 3 B° -> 7r°7r° events with a 
significance of 4.2 standard deviations including system- 
atic uncertainties. We measure a branching fraction 
B(B° -► 7r vr ) = (2.1 ± 0.6 ± 0.3) x 10~ 6 , where the 
first error is statistical and the second is systematic. The 
branching fraction is an average for B° and B° decays. 
The systematic uncertainties from PDF variations and 
efficiency have been combined in quadrature. This result 
is consistent with, and supersedes, our previous limit for 
this decay |l2| : it is also consistent with other prior lim- 
its 0| . The observed B° ~ * 7r°7r° branching fraction is 
larger than predicted theoretically. 




02468 10 02468 10 
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FIG. 1: The distributions of a) rriES, b) AE, and c) Fisher 
discriminant F for candidates in the signal data sample, and 
d) the F distribution for candidates in the sideband data sam- 
ple. The solid lines show the PDF for signal plus background. 
For ot.es and AE, the dashed lines show the PDF for the qq 
background. The abscissa in c) and d) is the F bin number, 
where the bins have been chosen so that each bin contains 
approximately 10% of the distribution for the B° — > D^nn 
data sample. 



TABLE I: A summary of systematic uncertainties listed as 
the change in the fitted event yield, Atijo^^.o, for different 
parameter variations. 



Parameter 



An f 



qq 7TtES shape parameter 
qq AE quadratic polynomial 

n B ±_ p ± n o 

B° -> ty°ty°AE mean 



±2.0 

+0.9 
-1.0 

±0.9 

+0.6 
-1.0 
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FIG. 2: The distributions of a) mES, b) AE, and c) Fisher 
discriminant F for candidates in the signal data sample that 
satisfy an optimized requirement on the signal probability, 
based on all variables except the one being plotted. The frac- 
tion of signal events included in the plots is 24%, 42% and 
74% for mES, AE, and F, respectively. The PDF projections 
are shown as a dashed line for qq background, a dotted line 
for B ± -> p ± -K°, and a dashed-dotted line for B° -> vrV 
signal. The solid line shows the sum of all PDF projections. 
The PDF projections are scaled by the expected fraction of 
events passing the probability-ratio requirement. The ratio 
— 2 In {C/Cmax) is shown in d) where the dashed line corre- 
sponds to statistical errors only and the solid line corresponds 
to statistical and systematic errors, as applied for the calcu- 
lation of significance. 



